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ABSTRACT 



B.JMeneux 7 , P. Oesch 1 , R. Pello 6 , C. 
Y. Taniguchi 23 , M. Zamojski 5 



We have used the zCOSMOS-bright 10k sample to identify 3244 Spitzer/MIPS 24 fim -selected galax- 
ies with 0.06 < S , 24 Aim < 0.50 mJy and I AB < 22.5, over 1.5 deg 2 of the COSMOS field, and studied 
different spectral properties, depending on redshift. At 0.2 < z < 0.3, we found that different redden- 
ing laws of common use in the literature explain the dust extinction properties of ~ 80% of our infrared 
(IR) sources, within the error bars. For up to 16% of objects, instead, the HaA6563 /H/3 A4861 ratios 
are too high for their IR/UV attenuations, which is probably a consequence of inhomogenous dust dis- 
tributions. In only a few of our galaxies at 0.2 < z < 0.3 the IR emission could be mainly produced by 
dust heated by old rather than young stars. Besides, the line ratios of ~ 22% of our galaxies suggest 
that they might be star- formation/nuclear- activity composite systems. At 0.5 < z < 0.7, we esti- 
mated galaxy metallicities for 301 galaxies: at least 12% of them are securely below the upper-branch 
mass-metallicity trend, which is consistent with the local relation. Finally, we performed a combined 
analysis of the H^ equivalent-width versus D n (4000) diagram for 1722 faint and bright 24 galaxies 
at 0.6 < z < 1.0, spanning two decades in mid-IR luminosity. We found that, while secondary bursts of 
star formation are necessary to explain the position of the most luminous IR galaxies in that diagram, 
quiescent, exponentially-declining star formation histories can well reproduce the spectral properties 
of ~ 40% of the less luminous sources. Our results suggest a transition in the possible modes of star 
formation at total IR luminosities Ltir ~ (3 ± 2) x 10 11 L©. 

Subject headings: infrared: galaxies - galaxies: evolution - galaxies: abundances - galaxies: starburst 
- galaxies: active 
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1. INTRODUCTION 

Investigating the composition of the extragalactic IR 
background (Puget et al. 1996; Elbaz et al. 2002; Dole 
et al. 2006) at z < 1 is of key importance to understand 
galaxy evolution and the global decline of the star forma- 
tion rate density in the last 8 Gigayears (Gyr). At z > 1, 
the IR background is mainly dominated by luminous and 
ultra- luminous IR galaxies (LIRGs and ULIRGs, with to- 
tal IR [5 - 1000 /im] luminosities 10 11 < L T ir < 10 12 L Q 
and Ltir > 10 12 L Q , respectively). However, these pop- 
ulations are progressively less important at lower red- 
shifts, becoming rare by the present epoch. The most 
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typical IR sources at z < 0.7 are IR normal galaxies, i.e. 
those with total IR luminosities Ltir < 10 11 L Q . Their 
study is fundamental to clarify the production of the bulk 
of the IR background in the second half of cosmic time. 

Our knowledge of IR galaxies at low redshifts was first 
possible thanks to studies conducted with the Infrared 
Astronomical Satellite (IRAS) and the Infrared Space 
Observatory (ISO) (see e.g. Aussel et al. 1999; Gen- 
zel & Cesarsky 2000, for a review), which revealed that 
both star formation and nuclear activity were associated 
with the dust emission. More recently with the advent of 
the Spitzer Space Telescope (Werner et al. 2004) and the 
Akari Telescope (Matsuhara et al. 2006), a much deeper 
insight into the nature of the z < 1 IR galaxy population 
has been possible. We have now characterised their lu- 
minosity evolution (e.g. Le Floc'h et al. 2005; Caputi et 
al. 2007; Magnelli et al. 2009), stellar masses (e.g. Caputi 
et al. 2006a, b), and spectral energy distributions (SEDs; 
e.g. Rowan- Robinson et al. 2005; Rocca-Volmerange et 
al. 2007; Takagi et al. 2007; Bavouzet et al. 2008; Kartal- 
tepe et al. 2009; Symeonidis et al. 2009), among others 
(see Soifer et al. 2008, for a review). 

Spectroscopic surveys at different wavelengths can add 
very valuable information that help understand the na- 
ture of the z < 1 IR sources. In addition to provid- 
ing accurate redshift determinations, spectroscopic data 
allow us to characterise other aspects of IR galaxies, 
such as their extinction properties, metallicities and gen- 
eral chemical composition. The initial optical/near-IR 
spectroscopic studies of IR galaxies at z < 1 targeted 
IRAS and ISO sources (e.g. Kim et al. 1995; Veilleux et 
al. 1995, 1999; Franceschini et al. 2003; Flores et al. 2004; 
Arribas et al. 2008; Garci'a-Marin et al. 2009). In the 
Spitzer era, and in the light of larger spectroscopic cam- 
paigns, more systematic studies of IR galaxies and their 
environments have been possible (e.g. Choi et al. 2006; 
Papovich et al. 2006; Caputi et al. 2008 [C08 hereafter], 
2009; Desai et al. 2008). 

In particular, optical spectra are useful to disentangle 
the power source of the IR emission, which generally oc- 
curs when the galaxy dust re-processes the energy pro- 
vided by the ultra-violet (UV) photons of young stars 
or active nuclei. However, in some cases, the observed 
IR emission might not be associated with any on-going 
stellar or nuclear activity but be rather produced by 
dust heated by old stars (e.g. Gordon et al. 2000; Kong 
et al. 2004; Cortese et al. 2006; but also see Young et 
al. 2009), or even direct stellar emission. The importance 
of these mechanisms within the IR galaxy population is 
not yet well established. The identification of 'passive' 
IR sources is important to quantify the degree of con- 
tamination in the IR-derived cosmic star formation rate 
densities. Spectroscopic data can be very helpful to elu- 
cidate the presence of star formation in sources that are 
presumably dominated by old stellar populations. 

Spectral features also carry the imprints of the galaxy 
star formation history. Thus, for IR galaxies, optical 
spectra can be used to investigate at what stage of the 
galaxy history the IR phase is produced. Locally, the 
rare LIRGs and ULIRGs are known to be mainly the 
consequence of galaxy mergers that induce bursts of star 
formation and nuclear activity (e.g. Armus et al. 2007). 
Instead, more typical IR galaxies at z = are usually 



less disturbed, and appear to form their stars in a more 
regular way. There are several clues that suggest that, in 
the past, when LIRGs and ULIRGs were more common, 
both the burst-like and quiescent modes of star formation 
could also lead to the IR phase (e.g. Elbaz et al. 2007; 
Daddi et al. 2008). In spite of these different pieces of 
evidence, it is not yet clearly established whether there 
exists a transition in the possible modes of star formation 
between galaxies of different IR luminosities. One of the 
main aims of this work is to study such a transition in 
IR galaxies at 0.6 < z < 1.0. 

The Cosmic Evolution Survey (COSMOS; Scoville et 
al. 2007) has been designed to probe galaxy evolution, 
star formation and the effects of large-scale structure up 
to high redshifts. The COSMOS field is defined by its 
Hubble Space Telescope/ Advanced Camera for Surveys 
(EST /AGS) 2 deg 2 coverage (Koekemoer et al. 2007). 
There are multiple follow up observations carried out 
in the COSMOS field, ranging from X-rays to radio 
wavelengths. Among the photometric imaging observa- 
tions, COSMOS includes the full and homogeneous cov- 
erage of the field with the Spitzer Infrared Array Camera 
(Spitzer /IR AC; Fazio et al. 2004) and the Multiband Pho- 
tometer for Spitzer (MIPS; Rieke et al. 2004), as part 
of Spitzer cycle-2 and 3 Legacy Programs (Sanders et 
al. 2007). 

COSMOS also comprises a large spectroscopic follow- 
up program named zCOSMOS (Lilly et al. 2007), that 
is being performed with the Visible Multiobject Spec- 
trograph (VIMOS; Le Fevre et al. 2003) on the Very 
Large Telescope (VLT). The zCOSMOS survey has two 
parts, one of which consists in the observation of an / < 
22.5 AB mag limited sample of 20,000 galaxies over 1.7 
deg 2 of the COSMOS field (zCOSMOS-bright). Here, we 
make use of the first half of spectra observed and anal- 
ysed in zCOSMOS-bright, which constitute the so-called 
'zCOSMOS-bright 10k sample' (Lilly et al. 2009). 

This paper presents the second part of a study we have 
carried out to analyse the optical spectral properties of 
24 /im -selected galaxies in the COSMOS field, using the 
zCOSMOS-bright 10k sample. In this second part, we 
analyse 3244 IR sources with 0.06 < £24 /un ^ 0.50 mJy, 
which constitute the largest optical spectroscopic sam- 
ple of mid-IR-selected galaxies analysed to date in the 
redshift range < z < 1. The results presented in this 
paper complement our first study of the optical spec- 
tra of brighter MIPS 24 /im sources (C08), and focus on 
the most outstanding aspects discussed in this previous 
work. The study of these new fainter galaxies is of par- 
ticular importance: they largely outnumber our previous 
sample and, especially, span the typical IR luminosities 
that dominate the IR background at z < 1 (Le Floc'h et 
al. 2005; Caputi et al. 2007; Magnelli et al. 2009), cover- 
ing a range in Ltir of e.g. a few times 10 10 to 9 x 10 10 L© 
at z = 0.5, and between Ltir ~ 9 x 10 10 and 5 x 10 11 L Q 
at z = 1. 

The layout of this paper is as follows: in Section Sj2j 
we describe the deep MIPS 24/xmand zCOSMOS-bright 
datasets and the cross-correlation of the two catalogues. 
In Section we give details about the rest-frame mid-IR 
luminosities for our galaxies. We present our results in 
Sections 2] and based on a variety of line diagnostics 
in different redshift ranges: we extensively discuss the 
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dust extinction properties of our galaxies at 0.2 < z < 
0.3, and study metallicities and stellar masses at 0.5 < 
z < 0.7. Later on, we perform a combined analysis of 
our present faint and previous brighter IR samples at 
0.6 < z < 1.0 to investigate the possible mechanisms of 
star formation in IR galaxies, and how they vary with 
IR luminosity. Finally, we summarise and discuss our 
results in Section £j6j We adopt throughout a cosmology 
with H = 70 km s^Mpc" 1 , fi M = 0.3 and Q A = 0.7. All 
quoted magnitudes and colours refer to the AB system. 

2. DATASETS 
2.1. The zCOSMOS-bright 10k sample 

The zCOSMOS survey (Lilly et al. 2007) is a very large 
spectroscopic program being performed with VIMOS (Le 
Fevre et al. 2003) on the VLT. VIMOS is a multi-slit 
spectrograph with four non-contiguous quadrants that 
cover 7x8 arcmin 2 each. zCOSMOS is designed to have 
a uniform pattern of pointings and a multiple-pass strat- 
egy, which guarantee a uniform coverage of the field. The 
technical characteristics of this survey are explained in 
further detail by Lilly et al. (2007). 

In zCOSMOS-bright, the vast majority of targets are 
selected randomly from a complete I < 22.5 AB mag cat- 
alogue (the 'parent catalogue'). This selection /-band 
magnitude corresponds to a total EST /AGS F814W 
magnitude in most cases. zCOSMOS-bright also includes 
a set of targets that have been selected for specific rea- 
sons, e.g. X-ray or radio sources. All the observations are 
performed with the R ~ 600 VIMOS medium-resolution 
(MR) grism, with a spectral coverage of (5500 — 9500) A 
and a dispersion of 2.55 A. The spectra are automat- 
ically reduced using the VIPGI software (Scodeggio et 
al. 2005), and the redshift determination of each source 
is manually checked on an individual basis by two inde- 
pendent reducers. 

The confidence in the redshift determination of each 
source is qualified with a flag, whose values can be: 4 
(completely secure redshift); 3 (very secure redshift, but 
with a very marginal possibility of error); 2 (a likely 
redshift, but with a significant possibility of error); 1 
(possible redshift); (no redshift determination); and 9 
(redshift based on a single secure narrow line, which is 
usually [Oil] A3727or HaA6563). The addition of +10 
to the flag indicates a broad-line active galactic nucleus 
(BLAGN). The confidence flags are manually assigned 
and reconciled by the two independent reducers (see Lilly 
et al. 2007, 2009 for further details). 

Independently, photometric redshifts for the entire 
zCOSMOS sample have been obtained with The Zurich 
Extragalactic Bayesian Redshift Analyzer (ZEBRA; 
Feldmann et al. 2006), using the UV through near-IR 
waveband data available for the COSMOS field. Also, a 
newer set of photometric redshifts based on 30-band pho- 
tometry have been derived by Ilbert et al. (2009) and Sal- 
vato et al. (2009). Tests performed on objects with dupli- 
cate spectra of different quality show that good-quality 
photometric redshifts as those obtained with the COS- 
MOS datasets can be useful to confirm less-secure spec- 
troscopic redshifts, such as those with flag=2, 1 or 9 (see 
Lilly et al. 2007). The differences between the photomet- 
ric and spectroscopic redshifts are Az < 0.08 X (1 + z) 
for > 93% of the objects with spectroscopic flag=2, 3 or 



4, and for 72% of those with a flag=l (Lilly et al. 2009). 

In this paper, we make use of the 'zCOSMOS-bright 
10k sample' over 1.5 deg 2 of the COSMOS field, in which 
10,580 out of 10,644 sources have / < 22.5 AB mag. 
This first set of zCOSMOS-bright spectra are now pub- 
licly available 24 . A full description of this dataset can be 
found in Lilly et al. (2009). 

2.2. The SCOSMOS MIPS deep 24 /mr catalogue 

The SCOSMOS survey (Sanders et al. 2007) comprises 
IRAC 3.6, 4.5, 5.8 and 8.0 /im and MIPS 24, 70 and 
160 /xm observations including the entire 2 deg 2 of the 
COSMOS- ACS field, and are part of the Spitzer cycle-2 
and 3 Legacy Programs. 

In cycle-3, the COSMOS field has been mapped 
at 24 /im down to a flux density S24 ^ m w 0.06 mJy 
(SCOSMOS-deep). The source extraction and pho- 
tometric measurements on the 24 /im maps have been 
performed with the SExtractor software (Bcrtin & 
Arnouts 1996) and the DAOPHOT package (Stet- 
son 1987), respectively. A point spread function (PSF) 
fitting technique for measuring the photometry has been 
necessary to deal with blending problems. Further de- 
tails on the SCOSMOS survey as well as the data re- 
duction and source detection procedures are given by 
Sanders et al. (2007) and Le Floc'h et al. (2009). 

The final SCOSMOS-deep 24 /xm catalogue contains 
52,092 sources with S24 fan > 0.06 mJy over the entire 
COSMOS field, with -70% of them lying within the 
zCOSMOS-bright lOk-sample area (see Figured]). 

2.3. Cross-correlation of the catalogues 

We performed the cross-correlation of the MJP^-deep 
and zCOSMOS-bright 10k catalogues in the same way 
as we did it for the shallow MIPS 24 /im sample in 
C08, so we only summarise the procedure here, quoting 
the figures corresponding to the matching of the deep 
24 /im sample. 

Firstly, as the vast majority of our zCOSMOS-bright 
10k sources are randomly taken from an / < 22.5 AB 
mag complete catalogue, it is important to make clear 
what fraction of the S24 Mm > 0.06 mJy sources can be 
identified with such an optical magnitude cut. To deter- 
mine this fraction, we cross-correlated the SCOSMOS- 
deep 24 /xm catalogue with the zCOSMOS-bright 'parent 
catalogue'. We obtained that around 39% of SW^m > 
0.06 mJy sources are associated with an / < 22.5 AB 
mag counterpart and, therefore, our analysis can be con- 
sidered as representative of this fraction of 24 /im sources. 

The cross-correlation of the deep 24 /im catalogue with 
the zCOSMOS-bright 10k catalogue allowed us to iden- 
tify 4497 IR sources, using a matching radius of 2 arcsec. 
This figure is the result of simultaneously considering: 
the fraction of all 24 /im sources lying in the zCOSMOS- 
bright lOk-sample field (—70%); the fraction of sources 
selected with the I < 22.5 mag cut (39%), and the 
zCOSMOS-bright 10k sampling rate, i.e. ~ 1/3 of the 
/ < 22.5 sources. 

To assess the presence of false identifications, we 
performed a separate cross-correlation of the deep 
24 /im catalogue with a deep complete / < 25 AB 

24 http:/ /archive. eso.org/cms/eso-data/data-packages/zcosmos- 
data-release-dr2 / 
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mag catalogue available for the COSMOS field (Capak 
et al. 2007). We found that a small fraction of our 
24 /im sources (138 out of 4497, i.e. 3%) have an opti- 
cal I < 25 AB mag association that is closer than the 
zCOSMOS-bright 10k counterpart. We excluded these 
138 sources from our analyisis. 

Out of the remaining 4359 24 /zm /zCOSMOS identi- 
fications, 3865 sources (i.e. ~89%) have a good-quality 
redshift determination (as indicated by their zCOSMOS 
quality flags; see Lilly et al. 2007 and C08). 3838 of these 
sources are galaxies, while the other 27 are galactic stars. 

Finally, our sample of 3838 24 (im galaxies with good- 
quality zCOSMOS-bright redshifts and spectra contains 
594 objects that we have already analysed in C08. (Note 
that the very small discrepancy of 2% in the number of 
bright 24 /zm galaxies analysed in C08 and found here is 
due to the update of the 24/zmand zCOSMOS-bright 
10k catalogues). In this work, we study new 3244 galax- 
ies with 0.06 mJy < S24 pm ^ 0.50 mJy, and combine 
or compare with our previous results for the brighter 
24 /im sample, where appropriate. 

Nearly 97 (99)% of the 3244 galaxies in our final sample 
lie at z < 1.2 (z < 1.5), as it can be seen in Figure[2j This 
figure also indicates the redshift bins and subset of galax- 
ies that we effectively consider for spectroscopic analysis 
in the following sections. The 1% high-redshift tail at 
z > 1.5 is almost exclusively composed of sources identi- 
fied as BLAGN with the zCOSMOS spectra. Within our 
24 izm /zCOSMOS sample, there are 88 BLAGN in total 
at any redshift. Roughly half of them are at z > 1.5. 

3. IR LUMINOSITIES 

We computed the rest- frame 24 /zm luminosities 
vL 2A Mm of the 3244 galaxies in our final sample in the 
same way as in C08: we used the observed 24 /zm fluxes, 
the zCOSMOS redshifts, and k-corrections derived us- 
ing the Lagache et al. (2003, 2004) IR spectral en- 
ergy distribution (SED) models. Our galaxies have 
24 /zm luminosities vL 2i ^ m w 10 s L Q at z — 0.1, roughly 
between 10 9 and 10 10 Lq at z = 0.5, and between 10 10 
and ~ 10 11 Lq at z = 1. 

The estimated total IR luminosities Ltir of these 
galaxies depend on the recipe used to convert monochro- 
matic into bolometric IR luminosities. By adopting 
the vL^^-Lttr calibration of Bavouzet et al. (2008), 
we obtain that our galaxies span Ltir values between 
a few times 10 10 and 9 x 10 10 Lq at z — 0.5, and 
L Tm w 9 x 10 10 to 5 x 10 n L Q at z = 1. These to- 
tal IR luminosities imply that our sample is composed of 
IR normal galaxies at intermediate redshifts, and mainly 
LIRGs by z = 1 . The analysis of IR sources in these lumi- 
nosity ranges is particularly important, as they make the 
bulk of the IR background at these redshifts (Le Floc'h 
et al. 2005; Caputi et al. 2007; Magnelli et al. 2009). 

In general, the SED models and recipes used to esti- 
mate IR luminosities apply only to IR galaxies dominated 
by star formation. The SED of some AGN, instead, are 
characterised by a power-law shape /„ cx v a (with a < 0; 
see e.g. Elvis et al. 1994; Alonso-Herrero et al. 2006). As 
in C08, we analysed the Spitzer/IRAC photometry avail- 
able for our galaxies to determine those cases in which 
the fluxes in the four IRA C channels 3.6, 4.5, 5.8 and 8.0 
/zm are consistent with a single power law, within the 
error bars. We found that 96 of our 3244 galaxies have 



an IRAC power-law SED, and nearly 80% of them are 
also zCOSMOS-classified BLAGN. 

4. SPECTRAL DIAGNOSTICS AT 0.2 < z < 0.3 
4.1. Balmer lines and optical colours 

The wavelength coverage (5500-9500) A of the VIMOS 
MR zCOSMOS spectra allows for different Balmer lines 
to be observed at different redshifts. The spectral qual- 
ity is quite uniform up to wavelengths A w 8000A, but is 
degraded by fringing at longer wavelengths, especially at 
A > 8500A. Thus, for example, although one can in prin- 
ciple observe the HaA65631ine up to redshift z w 0.45, 
reliable measurements of this line on individual spectra 
are generally only possible for sources at z < 0.3. 

We measured the fluxes and equivalent widths (EWs) 
of all lines in our spectra by direct integration on the 
rest-frame spectra. To minimise the errors introduced 
by fringing, we manually checked the measurements of 
those lines at observed 8000 < A < 8500 A which had 
particular problems, such as a significantly shifted cen- 
tre. We have not performed line measurements beyond 
observed A > 8500A. 

In this section, we focus on the analysis of those spec- 
tra in which the two main Balmer lines Ha A6563 and 
H/3A4861are simultaneously present, as their flux ra- 
tio provides information on the dust obscuration af- 
fecting the on-going star formation in IR galaxies. 
With this purpose, and willing to avoid measuring 
lines both in the fringing region and the blue edge of 
the spectra, we performed simultaneous Ha A6563 and 
H/3 A4861 measurements only for the 204 IR galaxies that 
lie at 0.2 < z < 0.3 within our sample (none of which is 
classified as a BLAGN). 

We corrected both the Ha A6563 and 
H/3 A4861 emission line measurements for stellar ab- 
sorption. We estimated the stellar absorption correction 
in each case by fitting the continuum of each spectrum 
with synthetic stellar SEDs from the Bruzual & Chariot 
2007 library (Bruzual & Chariot 2003; Bruzual 2007). 
We also applied aperture corrections to account for the 
fact that the slits in the VIMOS masks have a width 
of only 1 arcsec. To obtain the aperture corrections, 
we convolved each zCOSMOS spectrum in our sample 
with the Subaru Suprime-Cam R and /-band filters, 
and compared the resulting magnitudes with the total 
magnitudes of our galaxies, as available in the COSMOS 
optical photometry catalogues (Capak et al. 2007). 

Within our sample at 0.2 < z < 0.3, 189 out of 204 
(>92%) galaxies have Ha A6563 EW > 5 A (before stel- 
lar absorption correction) . This means that these galax- 
ies have an on-going star formation that is significant 
with respect to the underlying stellar population. The 
remaining 15 galaxies with Ha A6563 EW < 5 A will be 
analysed later, but we anticipate that only among them 
can we expect to have sources in which the IR emission 
is likely to be associated with dust heated by old stars 
rather than UV photons from on-going star formation. 
The adopted limit (EW = 5 A) to separate our galaxy 
sample is the same as the one we used in C08 and, of 
course, is arbitrary. However, no result or conclusion in 
this work depends on the exact value of that limit. 

Among the 189 24 ztm galaxies with HaA6563EW > 
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5 A at 0.2 < z < 0.3, 94 also have H/3A4861EW > 
5 A, while the other 95 have lower H/3A4861EW val- 
ues. This is not very surprising as, even in the absence 
of extinction, there is no reason for both Ha A6563 and 
H/3 A4861 EWs to be similar. Their ratio depends on 
the galaxy age and star formation history (because 
they determine the slope of the continuum between the 
HaA6563and H/3 A4861 wavelengths). In addition, any 
possible differential extinction between gas and stars 
could also influence the EW ratio. We discuss in more 
detail the extinction properties of these sources in Sec- 
tion OJ 

Although here we classify our galaxies according to 
their HaA6563and H/3A4861EWs (which measure the 
relative importance of the on-going and the past inte- 
grated star formation histories), it is instructive to see 
the behaviour of these quantities with respect to the line 
luminosities (Figure [3]). Globally, the EWs and luminosi- 
ties relate to each other, although with non-negligible 
dispersion. This is both due to the different optical con- 
tinuum levels in IR galaxies, and the spectral measure- 
ment errors. 

In this section, our main interest is to investigate 
whether IR galaxies with different B aimer line EWs can 
be recognised by their optical broad-band colours, and 
how, in general, the optical colour distribution of IR 
galaxies compares with the colour distribution of the en- 
tire zCOSMOS 10k sample. 

The left-hand panel of figure 2] shows the distribu- 
tion of the observed (B — i) colours for all zCOSMOS- 
bright 10k galaxies at 0.2 < z < 0.3, compared to 
the colour distribution for the 24 /im /zCOSMOS sources 
(note that the number of all zCOSMOS-bright galaxies 
is renormalised to 204 to coincide with the number of 
24 /im /zCOSMOS sources in this redshift range). On the 
right-hand panel, we plot the respective (B — i) colour 
distributions for the 24 fim /zCOSMOS sources with 
Ha A6563 and H/3 A4861 EW > 5 A; Ha A6563 EW > 5 A 
and H/3A4861EW < 5 A; and both EW < 5 A. In all 
cases, the (B — i) are total (i.e. aperture-corrected) 
colours measured on the Subaru Suprime-Cam images 
for COSMOS (Taniguchi et al. 2007). 

The comparison of these different histograms clearly 
shows that, while the overall zCOSMOS galaxy popula- 
tion at 0.2 < z < 0.3 displays a bimodal colour distri- 
bution, this bimodality is not seen for the IR sources. 
This is basically due to the relatively higher fraction of 
'green valley' 1.3 < (B — i) < 1.7 galaxies within the 
24 /im sample, and lower fraction of blue (B — i) < 1.3 
galaxies, which suggests that the zCOSMOS sources de- 
tected at 24 /im are on average more obscured. 

When we analyse the {B—i) colour distributions for the 
IR galaxy samples with different values of Ha A6563 and 
H/3 A4861 EWs, we see the following: those galaxies with 
both HaA6563and H/3A4861EW > 5 A make most (~ 
75%) of the optically-blue [B — i) < 1.3 IR galaxies. In- 
stead, the 'green valley' and red IR galaxies are mainly 
sources with H/3A4861EW < 5 A (58% and 78%, re- 
spectively). Therefore, although there is no one-to-one 
relation between the line EWs and the broad-band flux 
ratios, one can roughly predict the spectral character- 
istics of IR galaxies according to their different optical 
colours. 



A similar conclusion can be extracted by plotting the 
observed (B — i) colours as a function of the Balmer 
decrement L(H a) / 'L(H /3) (Figure [5]). In this diagram, it 
is clear that galaxies with H/3 A4861 EW < 5 A dominate 
the region defined by (B — i) > 1.3. The vast majority 
of our 24 /im galaxies display Ha A6563 /H/3 A4861 values 
above the intrinsic ratio in the typical case B recombina- 
tion with T = 10, 000 K (i.e. Ha A6563 /H/3 A4861 =2.87; 
Osterbrock 1989), independently of their optical colours, 
indicating in all cases significant dust attenuations. How- 
ever, nearly all the sources with the highest dust extinc- 
tions HaA6563/H/3A4861 > 7 have (B - i) > 1.3 (cf. 
Section 34~3)) . 

The stacked spectra of the sources with 
H/3A4861EW < 5 A, classified according to their 
(B — i) colours, are useful to summarise the average 
spectral properties (Figure These stacks show 

that the Balmer line EWs become smaller for redder 
colours, with the (uncorrected) H/3A4861EW ranging 
from (3.48 ± 0.50) A for the (B — i) < 1.3 stack, to 
(0.62 ± 0.20) A for that with (B - i) > 1.7. The average 
stellar-absorption-corrected L(H a) / 'L(H /3) ratio varies 
from (4.12 ± 0.52) to (10.3 ± 2.6), respectively. 

4.2. IR galaxies with Ha EW < 5 A 

Our 24 /im /zCOSMOS galaxy sample at 0.2 < z < 
0.3 contains 15 sources whose individual Ha A6563 EWs 
are smaller than 5 A. This indicates that either the on- 
going star formation in these galaxies is insignificant with 
respect to their integrated past star formation history, or 
that there is a very different degree of extinction in the 
gas associated with the new star formation and the old 
stars responsible for the continuum light. In this section, 
we investigate in detail the properties of these galaxies, 
and discuss in particular the origin of their IR emission. 

Figure [7] shows the average stacked zCOSMOS spec- 
trum of these 15 galaxies. This spectrum is charac- 
terised by a weak Ha A6563 line in emission with EW = 
(4.8 ± 1.0) A, and an H/3 A4861 line in absorption that in- 
dicates the importance of the underlying old stellar pop- 
ulations. The presence of prominent absorption Mgl line 
at A = 5175.4 A and NaD line at A = 5892.5 A confirms 
the dominance of the old stellar component. 

To better understand the nature of our 
Ha A6563EW < 5 A sources, we inspect their EST /ACS 
F814W-band morphologies on an individual basis. The 
AGS images reveal a variety of cases: some of these 
galaxies are irregular with clumps, favouring the case 
of heavily obscured star formation. A few are simply 
edge-on galaxies, so the large extinction is mainly an 
orientation effect. 

Five out of 15 of these sources are ellipticals or lenticu- 
lar galaxies. In principle, there is no indication of signif- 
icant star formation in these objects except for their IR 
emission. Only two of them are detected in the near-UV 
with NUV < 25.8 mag (cf. Section O . These galaxies 
are probably the only candidates in our 0.2 < z < 0.3 
sample for which the IR emission could be mainly asso- 
ciated with dust heated by the UV photons of evolved 
stars. 

Figure [8] shows the rest-frame broad-band SEDs of 
these five galaxies with presumably no star formation. 
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For reference, we added the SED model templates of 
three galaxy types: Elliptical, SO and Sc (from Pol- 
letta et al. 2007; based on the GRASIL code, Silva et 
al. 1998). For the spiral models, the templates in the 
wavelength range 5-12 /mi have been constructed using 
empirical spectra observed with ISO and Spitzer (see Pol- 
lctta et al. 2007 for more details). 

The comparison of our galaxy SEDs with the different 
models shows that, in all cases, the mid-IR emission is in 
excess of what is expected from pure stellar emission (as 
in the case of an elliptical galaxy model) . This indicates 
that the presence of dust is indeed necessary to explain 
the mid-IR emission in our galaxies. Note, however, that 
in the two middle-panel SEDs the mid-IR excess is small, 
i.e. the direct stellar emission constitutes a significant 
part of the 24 /xm flux. 

The IR emission of none of these five sources seems to 
be powered by an AGN, as none of them has an IRAC 
power-law SED or is detected in the X-ray catalogues 
available for the COSMOS fields (Hasinger et al. 2007; 
Elvis et al. 2009). The mid-IR luminosity distribution 
of these galaxies is similar to that of IR sources with 
larger-EW Balmer lines at similar redshifts. 

4.3. Dust extinction 

As in C08, we analyse our 24 /im galaxy sam- 
ple with simultaneous HaA6563and H/3A4861 
measurements at 0.2 < z < 0.3 to investigate the 
relation between the Balmer decrement L(Ha)/L(H(3) 
and the far-IR-to-ultra-violet (UV) luminosity ratio. 
In general, for star forming galaxies, the far-IR/UV 
ratio gives a measurement of the dust extinction ^2000: 
(£70 M m + L 200 o)/L 20QO « i 0-4A 2 ooo (e.g. Buat & 
Xu 1996; Gordon et al. 2000; Buat et al. 2002). 

Our results for bright IR galaxies with Ha A6563 and 
H/3 A4861 EW > 5 A at 0.2 < z < 0.3 (C08) showed that, 
for most of them, the Balmer decrement versus the far- 
IR/UV attenuation was compatible with the Calzetti et 
al. (2000) reddening law, within the error bars. Here, we 
would like to investigate whether this conclusion is still 
valid when considering fainter IR galaxies. 

We estimated the rest-frame 70 fim luminosity of our 
24 (im galaxies at 0.2 < z < 0.3 using their rest-frame 
24 /im luminosities and the vL 2 ^ Mm -vL„ A ' m conversion 
formula calibrated by Bavouzet et al. (2008). To ob- 
tain UV luminosities at rest-frame A = 2000 A, we used 
the GALEX near-UV data available for the COSMOS 
field (P.I.: D. Schiminovich), which has an observed effec- 
tive wavelength A e ff, = 2310 A. For NUV non-detected 
sources, we adopted a lower limit of 25.8 AB mag. 

Figure [9] shows the Balmer decrement L(Ha)/L(H(3) 
versus the extinction log 10 ([L 70A1 m + ^20oo]/^20oo), for 
our galaxies with HaA6563EW > 5 A at 0.2 < z < 0.3 
(circles and diamonds for H/3A4861EW > and < 5 A, 
respectively) . The values obtained for brighter IR galax- 
ies are also shown (crosses). Note that Figure [9] shows 
more bright 24 fim galaxies than those analysed in C08, 
as here we want to study all galaxies with Ha A6563 EW 
> 5 A (independently of the H/3 A4861 EW), as we do for 
our current fainter sample. This allows us to extend our 
dust extinction analysis to the most obscured sources. 

We see that, when considering all bright and faint 
24 /mi galaxies with HaA6563EW > 5 A, the spread in 



the extinction diagram is quite large. While many of our 
sources lie above or close to the relation derived from 
the Calzetti et al. (2000) reddening law, several others 
have lower Balmer decrements for their Lir/Ljjv atten- 
uations. The attenuation of some of these galaxies seems 
more compatible with the Small Magellanic Cloud red- 
dening law (SMC; Prevot et al. 1984), but unfortunately 
it is difficult in general to decide among the different red- 
dening laws given the line-ratio error bars. 

Around 5% of our 24 /mi galaxies are characterised by 
a L(Ha)/L(H0) ratio close to a T = 10, 000 K case B 
recombination, but have a substantial Lir/Luv attenu- 
ation. This effect has also been discussed in C08. These 
are usually cases of sources with an inhomogeneous dust 
distribution, where the VIMOS 1-arcsec slit only receives 
the optical light of the practically unobscured central 
region, while the large («6-arcsec full-width half max- 
imum) Spitzer /MIPS point spread function collects the 
IR light of the dustier surroundings. However, in our 
new fainter sample, we identify a more interesting prop- 
erty for several of these outlier galaxies: their line ratios 
suggest that they could be composite star-forming/ AGN 
systems (see Section §4.4|) . For these sources, trying to 
understand their dust attenuation with reddening laws 
calibrated on pure star-forming galaxies is probably mis- 
leading. 

On the other hand, we see that other galaxies in our 
sample are characterised by very large Balmer decre- 
ments: around 16% of them are at least la away from any 
of the common reddening laws (this figure includes ~13 
and 19% of sources with H/3 A4861 EW > 5 A and < 5 A, 
respectively). The use of any of these reddening laws 
with the corresponding Lir /Luv attenuations would sig- 
nificantly under-estimate the H/3 A4861 extinctions ob- 
served in these galaxies. Interestingly, very few of 
them have line ratios characteristic of composite systems 
( q4.4[) , so the plausible presence of obscured nuclear ac- 
tivity does not appear to be related to the high Balmer 
decrements. 

Among the brighter IR galaxies in the C08 sample, and 
taking into account all the sources with Ha A6563 EW > 
5 A, we find that around ~11% have Balmer decrements 
that are significantly above the predictions of the various 
reddening laws. 

A possible caveat is that these outliers could be sim- 
ply produced by an under-estimation of the stellar ab- 
sorption corrections applied to the H/3 A4861 lines, whose 
median is (3.1 ± 0.3) A and maximum value is 3.7 A. 
Such a correction is comparable to the same un-corrected 
H/3A4861EWs in some cases. To test this possibility, we 
calculated the hypothetical stellar absorption corrections 
that would be needed in order to make the outliers have 
Balmer ratios within lcr of the Calzetti et al. (2000) red- 
dening law. For this to happen, we found that around a 
half of the outliers would need stellar absorption correc- 
tions of more than 3.7 A, the maximum value obtained 
for any of our galaxies by fitting their spectral contin- 
uum. This fact suggests that, at least for half of them, 
the large displayed Balmer decrements are not an arti- 
fact of under-estimated stellar absorptions. Note that, 
in this exercise, we have not considered any variation 
in the Ha A6563 stellar absorption (in reality, the two 
corrections would increase simultaneously). Thus, our 
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test determines the maximum possible fraction of out- 
liers that could be corrected just by applying a larger, 
though still realistic, stellar absorption correction. 

A likely explanation for the highly extincted galax- 
ies is that they contain heavily obscured starbursts 
which are concentrated in their central regions. While 
(£70 + ^20oo)/-^20oo is related to the total amount 
of galaxy dust, L{Ha)/L{Hf}) is more sensitive to the 
dust concentration. So, even if the total amount of dust 
is not very large, the HaA6563 /H/3 A4861 ratio will be, 
when the dust is too concentrated towards the centre (the 
HaA6563and H/3 A4861 fluxes correspond to the central 
1-arcsec width light collected in the VIMOS slit) . Note 
that basically all of these galaxies are detected in the 
NUV, so either they are intrinsically bright at these wave- 
lengths or, indeed, the highly obscured region is very 
compact and the UV photons can leak from the less ob- 
scured surrounding regions. 

It is important to remark that the discrepancies be- 
tween the extinctions derived from the Balmer decre- 
ments and the IR/UV are expected to be less important 
at higher redshifts, because the spectral slits collect a 
more representative part of the total galaxy light. 

4.4. The BPT diagram 

The [OIII] A5007 /H/3 A486 1 versus 

[Nil] A6584 /Ha A6563 diagram (Baldwin, Phillips & 
Terlevich 1981; BPT hereafter) provides a useful di- 
agnostics to assess the diversity of the nature of a 
galaxy population. For example, one can disentan- 
gle the existence of nuclear activity from the high 
[Nil] A6584 /Ha A6563 ratios, or select candidates to 
star-forming/ AGN composite systems (e.g. Kauffmann 
et al. 2003a). 

Figure [TO] shows the BPT diagram for our 
24 /im galaxies at 0.2 < z < 0.3. Symbols are the 
same as in Figure [9] circles indicate the line ratios 
of galaxies with both HaA6563and H/3A4861EW 
> 5 A, while diamonds correspond to galaxies with 
HaA6563(H/3A4861) EW > (<)5A and different 
(B — i) colours. 

In all cases, the Ha A6563 and H/3 A4861 values are cor- 
rected for stellar absorption. However, we do not in- 
clude extinction corrections. Given the dispersion of 
these sources in the dust attenuation diagram presented 
in Figure [5J we preferred to present the crude line ra- 
tios in the BPT diagram. The errors of this approach 
are basically negligible ( < 0.0025 and < 0.02 on log 10 
of [Nil] A6584 /Ha A6563 and [OIII] A5007 /H/3 A4861 , re- 
spectively). 

Inspection of Figure [TD] shows that 24 ^m sources oc- 
cupy quite different places within the BPT diagram, pro- 
viding evidence of the heterogeneous character of the IR 
normal galaxy population. We obtained a similar con- 
clusion in C08 for brighter IR sources. 

In our present faint IR sample, we find a significant 
fraction of galaxies lying in the composite region of the 
BPT diagram. This fact suggests that simultaneous star 
formation and nuclear activity could be present in these 
sources. These galaxies make around 22% of our sources 
with HaA6563EW > 5 A at 0.2 < z < 0.3, including - 
24% of the sample with H/3 A4861 EW > 5 A, and 20% of 
the sources with H/3A4861EW < 5 A. In addition, one 



of the red galaxies with H/3 A4861 EW < 5 A is identified 
as a Seyfert galaxy. 

Within the brighter IR sample (crosses in Figure [TO]), 
the percentage of composite sources is < 15% (and quite 
smaller among the sources with H/3A4861EW > 5 A 
analysed in C08), and no Seyfert galaxy has been found. 
Among all the optically-selected zCOSMOS-lOk galaxies 
at similar redshifts, the percentage of composite galaxies 
in the BPT diagram is only ^10-12%. 

A bit more than a half of the composite galaxies in our 
faint IR sample are optically blue ((B — i) < 1.3), while 
the rest have redder colours. We note that, in spite of 
being spectroscopically classified as composites, none of 
these galaxies has an IRAC power-law SED or is detected 
in X-rays. This fact suggests that the plausible AGN 
component is not dominant, and that aperture effects 
could also have a role in the fraction of sources identified 
as composite. 

Indeed, if the AGN component were concentrated to- 
wards the centre, the line ratios observed in the 1 arcsec- 
wide optical spectra would not be the same as those 
obtained from the total galaxy light. In this case, one 
should probably expect to identify a lower fraction of 
composite systems among similar IR sources at higher 
redshifts, for which the central light would not be spec- 
troscopically resolved. 

5. SPECTRAL DIAGNOSTICS AT z > 0.5 
5.1. The mass-metallicity relation 

Atz> 0.5, the Ha A6563 and [Nil] A6584 emission lines 
are out of the coverage of the zCOSMOS-bright VIMOS 
spectra, but the [Oil] A3727 line starts to be visible in- 
stead. In particular, between z = 0.5 and z = 0.7, all the 
three [Oil] A3727 , H/3A4861and [OIII] A5007 emission 
lines are located in good regions of the zCOSMOS spec- 
tra (i.e. outside the regions severely affected by fring- 
ing). The flux line ratios of these three lines can be 
used to estimate the galaxy oxygen abundances (with e.g. 
the so-called R23 parameter, i.e. R23 = ([Oil] A3727 + 
[OIII] A4959, 5007) /H/3; see Pagel et al. 1979). 

The relation between galaxy mass and metallicity at 
different redshifts has been widely investigated from the 
observational (e.g. Lequeux et al. 1979; Zaritsky et 
al. 1994; Tremonti et al. 2004; Maiolino et al. 2008) and 
theoretical (e.g. de Lucia et al. 2004; Tissera et al. 2005; 
de Rossi et al. 2007; Dave & Oppenheimer 2007; Finla- 
tor & Dave 2008) points of view. The basis of this rela- 
tion is that galaxy growth and chemical enrichment are 
closely associated. However, during an episode of star 
formation, such as those characterising most IR galax- 
ies, gas inflows and outflows can alter the composition 
of the inter-stellar medium. This fact contributes to the 
scatter observed in the mass-metallicity relation and, in 
particular, to the existence of outliers. 

In C08, we found that around 1/3 of the LIRGs and 
ULIRGs for which we had metallicity measurements at 
0.5 < z < 0.7 were securely below the mass-metallicity 
relation. Our aim here is to investigate the fraction of 
outliers in the mass-metallicity relation among less lumi- 
nous IR sources at similar redshifts. Our present faint 
24 ^m sample at 0.5 < z < 0.7 spans around a decade in 
IR luminosities 3 x 10 10 < L TIR < 3 x 10 n L Q , i.e. it 
is mainly composed of IR normal galaxies and the less 
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luminous objects among the LIRGs. 

We derived metallicities for our 301 IR galaxies 
with H/3 A4861 EW > 5 A and positive [Oil] A3727 and 
[OIII] A5007 fluxes. We used the algorithm described by 
Maier et al. (2005), which is based on the Kewley & Do- 
pita (2002) model. We applied this algorithm on com- 
pletely corrected emission lines, i.e. corrected for stellar 
absorption, aperture effects and extinction. zCOSMOS 
BLAGN have been excluded from the analysis. 

To correct our lines for extinction, we assumed that 
the sum of the star formation rates (SFRs) derived from 
the UV and IR should be equal to that derived from 
the H/3 A4861 luminosity, as obtained applying formulae 
commonly used in the literature (Kennicutt 1998; see 
C08 for more details). This procedure could be dubious 
for some galaxies at low z < 0.3 redshifts: we remind 
that we have concluded in Section N4.3I that the Balmer 
decrements and IR/UV attenuations of up to 16% of our 
galaxies at 0.2 < z < 0.3 could not be reconciled with 
any of the commonly used reddening laws. 

However, the galaxies in our present sample are at 
somewhat higher redshifts (0.5 < z < 0.7), where the 
aperture effects should be less important. Based on this, 
and the fact that the typical reddening laws do roughly 
predict the extinctions of most IR galaxies even at low 
z, we decided to apply the SFR balance procedure - 
SFR(UV+IR)«SFR(H /3 )- to estimate the extinction of 
our IR galaxies at 0.5 < z < 0.7. In any case, as we 
discuss later, this assumption does not affect any of our 
conclusions. 

Figure [TT1 shows the [OIII] /[Oil] ratios versus the R23- 
parameter values for the 301 IR galaxies in our metallic- 
ity sample. The resulting metallicity versus stellar-mass 
diagram is shown in Figure 1121 Stellar masses for all our 
galaxies have been derived through the synthetic tem- 
plate fitting of their broad-band SEDs, from UV through 
near-IR wavelengths (Bolzonella et al. 2009), and cor- 
respond to a Salpeter (1955) initial mass funtion over 
stellar masses M = (0.1 - 100) M . 

Around 79% of galaxies in our sample have degen- 
erate O/H values (small blue circles and crosses in 
Figure 1121 showing the upper and lower alternative 
metallicities, respectively). The metallicity is in gen- 
eral two-valued with respect to the oxygen abundances 
[Oil] A3727 /H/3 A4861 and [OIII] A5007 /H/3 A4861 , be- 
cause low line ratios can be produced either by low abun- 
dances, or high abundances with very efficient line cool- 
ing. 

The remaining ^21% of our sample have non- 
degenerate metallicity solutions (large red circles in Fig- 
ure [12]). Not surprisingly, these galaxies are those with 
the largest and smallest R23 values (see Figure [11]) , as 
the oxygen-abundance versus metallicity curve is single- 
valued for extreme R23 (this is particularly true for 
the Kewley & Dopita calibration; see Kewley & Do- 
pita 2002). 

From the analysis of Figure 1121 we can extract two 
main conclusions. Firstly, the upper-branch metallicities 
of our sources with degenerate solutions -whose median 
redshift is z w 0.62- concentrate around the local mass- 
metallicity relation. As a reference for the local mass- 
metallicity relation, we considered the results obtained 
on Sloan Digital Sky Survey galaxies at z < 0.1 (Kewley 



& Ellison 2008; solid line in Figure [£2]). The metallicities 
in this local relation also correspond to the Kewley & 
Dopita (2002) calibration. We note that we could not 
have concluded on this effect in C08, due to the small 
number of sources with metallicity measurements in this 
previous work. 

Secondly, more than a half of our sources with non- 
degenerate metallicity solutions lie well below the lo- 
cal mass-metallicity relation, around which the upper- 
branch solution of our sources with degenerate metal- 
licities are located. Most of the non-degenerate metal- 
licities concentrate around the turning point of the 
oxygen-abundance versus metallicity R23 curve, i.e. at 
12 + log 10 (O/H) w 8.6, consistently with them having 
large R23 values. 

The sources with low non-degenerate metallicities dis- 
play a variety of stellar masses, putting a secure lower 
limit of ~ 12% for the outliers to the mass-metallicity 
relation within our sample at 0.5 < z < 0.7. This per- 
centage is considerably smaller than that estimated for 
more luminous IR galaxies at similar redshifts (C08), but 
we emphasize that it is only a lower limit, as there could 
be more low-metallicity sources among those with degen- 
erate abundances. The under-abundance is probably not 
the consequence of a line ratio produced by the presence 
of an AGN within our sources, as only two of these are 
detected in X-rays (plus- like symbols in Figure [T2]) . 

We observe that a significant outlier ratio to the mass- 
metallicity relation has also been reported in optically- 
selected samples at z < 1, as e.g. the Canada France 
Redshift Survey (Lilly et al. 1995) galaxies analysed by 
Maier et al. (2005; triangles in Figure [T2J see also Lilly 
et al. 2003). The stellar masses of Maier et al. sources 
are comparable to those of our 24 /im galaxies, suggesting 
that the ~ 12% outlier fraction we derived here is very 
likely a lower limit. 

To test whether the main conclusions obtained in 
this Section depend on the adopted dust-extinction- 
correction recipe, we derived new metallicity values for 
our galaxies based on the emission lines un-corrected 
for dust extinction, and compared with our previous re- 
sults. We found that the percentage of sources with non- 
degenerate metallicity values in the un-corrected case is 
19%, and the upper branch solutions in the degenerate 
cases are still compatible with the local mass-metallicity 
relation within the error bars. The percentage of secure 
outliers to this relation is ~ 9%. These figures are very 
similar to those obtained using the extinction-corrected 
emission lines. 

5.2. The IR phase in the galaxy star-formation history 
5.2.1. Individual D n (4000) measurements 

The 4000 A-break in galaxy spectra is produced by the 
differential contribution of young and old stars that, re- 
spectively, dominate the spectral light bluewards and 
redwards of A = 4000 A. The strength of this feature 
D n (4000), usually defined as the ratio of the average flux 
densities between the (4000 - 4100) and (3850 - 3950) A 
bands (Balogh et al. 1999), is a very good indicator of 
the time passed since the last significant burst of star 
formation occurred in the galaxy. 

Besides, the presence of higher-order Balmer lines in 
absorption, such as H5 A4102, indicates recent star for- 



The optical spectra of 24 /im galaxies in COSMOS: II 



9 



mation activity, which typically has occurred within the 
last Gyr (e.g. Dressier & Gunn 1983). The combined 
diagnostics of the D n (4000) parameter and Us EW in 
optical spectra is, then, commonly used as a powerful 
tool to constrain the galaxy star formation history (see 
e.g. Kauffmann et al. 2003b; C08). 

Our results in C08 showed that the vast majority of 
optically bright ULIRGs and some of the most luminous 
LIRGs at 0.6 < z < 1.0 were characterised by spec- 
tra with D n (4000) < 1.2. We interpreted this effect as 
the consequence of secondary bursts of star formation 
produced in galaxies with underlying old stellar popula- 
tions. Here, we do a combined analysis of the IR-brightcr 
C08 sample and our current fainter sample to investigate 
whether this mode of star formation is similarly valid for 
galaxies of different IR luminosities. 

Our present 24 /im sample at 0.6 < z < 1.0 is mainly 
composed of galaxies with rest- frame 24 /jm luminosities 
5 x 10 9 < vL 2 u A » m < 10 n L o , i.e. total IR luminosities 
5.3 x 10 10 < Ltir < 5 x 10 11 L Q . This range is partic- 
ularly important for IR galaxies at 0.6 < z < 1.0, as it 
contains most of the IR luminosity density at these red- 
shifts (the total IR luminosity function is characterised 
by L* w 2.5 x 10 11 L at z = 1; see Caputi et al. 2007). 
The brighter C08 sample extends the IR luminosity range 
to uLl 4 ^ ks 5 x 10 11 Lq, i.e. L Tm « 1.5 x 1O 12 L . 
Thus, the combined analysis we perform here covers two 
decades in mid-IR luminosities (or 1.5 decades in -Ltir)- 

We measured the D n (4000) parameter values on the 
spectra of a total of 1722 IR sources at 0.6 < z < 1.0. 
This figure excludes BLAGN. The resulting vlA A ^ m ver- 
sus D n (4000) diagram is shown in Figure [131 In this 
figure we can see that, in contrast to the more luminous 
IR galaxies studied in C08, the less luminous sources 
display a wide range of D n (4000) values, from « 0.8 to 
w 1.8. Still, a bit more than 60% of them have spectra 
with D n (4000) < 1.2. 

We used the diagram in Figure [13] to classify our 
galaxies in different bins of D n (4000) values and 
z/L^Mm mmmosities We gtacked the zCOSMOS spec- 
tra of all galaxies in each bin, in order to measure their 
average spectral properties (cf. Section £15.2. 2\ . The bins 
have been chosen in such a way to produce a good sam- 
pling of the (D n (4000); vL 2 ^^ ) space with a sufficient 
number of objects. 

The D n (4000) values shown in Figure [TBI and used for 
the classification are not corrected for dust extinction. 
This assures that we stacked spectra with similar char- 
acteristics prior to any dust extinction assumption. In 
any case, the impact of the dust extinction corrections is 
small: for example, an extinction Ay = 2 only produces 
a correction of ~8% on the D n (4000) values. 

5.2.2. The H s EW - D„(4000) diagram: discerning the 
modes of star formation 

Locating our galaxies in the H 5 A4102 EW - D n (4000) 
diagram allows us to have a better insight of the time at 
which the IR phase occurs within the galaxy star forma- 
tion history. 

As we explained in C08, the Us absorption EWs in 
the zCOSMOS spectra of IR galaxies are usually quite 
small and the associated errors too large as to rely on 
the individual measurements. Therefore, to perform 



EW measurements with more precision, we stacked the 
IR galaxy spectra. 

We classified all (bright and faint) 24 /im galaxies at 
0.6 < z < 1.0 in 14 groups according to their D n (4000) 
and vL^ A ^ m values, and averaged their respective spec- 
tra. We then measured the D n (4000) parameter and H5 
EW on each composite spectrum. The number of galax- 
ies in each stack varies between 495 and 28. We have 
not separated our galaxies according to their redshifts 
because the binning in the 0.6 < z < 1.0 redshift range 
does not produce any relevant effect in all the analysis 
presented here. 

The resulting H 5 EW - D n (4000) diagram for our 
stacked spectra is shown in Figure [14J Positive EW val- 
ues in this diagram indicate in absorption, and are 
corrected for line filling. To do this correction, we mea- 
sured the H7 emission line on each composite, corrected 
it for stellar absorption, and assumed an intrinsic decre- 
ment ITy/H(5 = 1.80 (for a case B recombination with 
temperature T = 10, 000K; Osterbrock 1989). The final 
H,5 EW values carry significant error bars, mainly be- 
cause the crudely measured values are quite small and 
also because of the uncertainties in the corrections. The 
D n (4000) values shown in Figure [141 do include dust ex- 
tinction corrections, which we estimated using the me- 
dian extinction of the galaxies in each stack group. The 
total IR luminosities shown in the labels have been ob- 
tained by converting the limits of the vl?^ Mm bins with 
Ltir = 6858 x (vLl A ^ m f- 71 (Bavouzet et al. 2008). One 
should keep in mind that using a different recipe to con- 
vert mid-IR into total IR luminosities will produce some 
variations in the derived Ltir values. 

The tracks shown with thick solid lines in Fig- 
ure [14] correspond to galaxy models with different 
exponentially-declining star formation rates SFR oc 
e^ t ' T ', where r is the characteristic decaying time that 
varies between r = 0.01 Gyr and r — > 00 (i.e. a model 
with a constant SFR). We have produced these tracks 
by measuring the Us EW and D n (4000) parameter on 
galaxy synthetic SEDs generated with the GALAXEV 
code and the 2007 version of the Bruzual & Chariot mod- 
els (Bruzual k Chariot 2003; Bruzual 2007). All these 
tracks are plotted only up to ages < 6 Gyr, as this is 
approximately the age of the Universe at z = 1 . 

As we have discussed in C08, the average position of 
the most luminous IR galaxies in the EW - D n (4000) 
diagram cannot be reproduced by any of these simple r- 
decaying star formation histories. Our present analysis 
confirms that this is the case for galaxies with total IR 
luminosities Ltir > 3 x 10 11 L (squares with the two 
largest sizes in Figure [14]) . Instead, secondary bursts of 
star formation in galaxies with underlying old ( > 3 Gyr) 
stellar populations appear as a plausible mechanism to 
explain the spectral characteristics of these most lumi- 
nous IR sources (thin solid lines in Figure [14]). We note 
that bursts producing 5-10% of the galaxy stellar mass 
are sufficient to drastically move the position of a few- 
Gyr-old galaxy in the diagram from D n (4000) > 1.5 to 
D n (4000) < 1.2. 

Interestingly, although the average (H<5 EW;D n (4000)) 
points for the most luminous IR galaxies with 
D n (4000) < 1 do seem to be consistent with r-decaying 
models, they intersect their tracks at ages < 0.1 Gyr. 
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As it is quite unlikely that there exist many of such 
young galaxies at 0.6 < z < 1.0, the scenario of sec- 
ondary bursts of star formation is still the most reason- 
able to explain the properties of these sources. This ar- 
gument is supported by the fact that the stacked spectra 
of these galaxies show some signatures of old underly- 
ing stellar populations, such as the presence of Call H 
& K absorption lines. Also, these sources are known 
to have already assembled intermediate-to-large stellar 
masses (M > 1 x 10 10 M Q ; see e.g. Caputi et al. 2006a), 
suggesting that we are not witnessing their first episode 
of star formation. 

As we discussed before, for the less luminous IR 
(^tir < 3 x 10 11 L Q ) galaxies (squares of the two small- 
est sizes in Figure IH)). the D n (4000) parameter has a 
wider range of values. Taking into account their average 
position in the Bbj EW - D n (4000) diagram and other 
features in the stacked spectra, we can obtain clues on 
the star formation histories of the sources in different 
D n (4000) bins. 

The sources with D n (4000) < 1 are characterised by 
an average spectra with relatively weak absorption lines, 
and an Us line with EW < 2 A (that, in the case of 
the less luminous IR galaxies, is actually observed in 
emission; see upper panels in Figure [T5|) . They also 
show a prominent [Oil] emission line and, those with 
^tir < 1-7 x 10 11 L Q , a modest but significant [Nelll] 
line. These spectral characteristics, and the resulting av- 
erage position of these sources in the Ebj EW -D n (4000) 
diagram, suggest that these are either really young galax- 
ies or, more likely, they are experiencing a major sec- 
ondary burst of star formation that hides any plausible 
older stellar component. Note that the ionisation lines 
could also be the consequence of the presence of AGN. 
However, only one of these sources is detected in the 
X-ray catalogues for the COSMOS field, and only two 
have an IRAC power-law SED (one of which is the X-ray 
source). 

The sources with Ltir < 3 x 10 n L Q in the 1.0 < 
D n (4000) < 1.2 bin are characterised by quite differ- 
ent average stacked spectra (lower panels in Figure [15|) . 
These galaxies display prominent Call H & K absorption 
lines, indicating that the underlying old stellar compo- 
nent is important. Their average position in the diagram 
at (extinction-corrected) D n (4000) ks 1.02 seems compat- 
ible with some of the T-decaying models at young (< 1 
Gyr) ages. However, the optical spectra indicate that 
these galaxies are very probably older. Thus, once more, 
secondary bursts are the most plausible mechanism for 
the on-going star formation in these objects. 

The stacking points indicating Ltir < 3 x 10 11 L Q 
galaxies at D n (4000) > 1.15 in Figure fT4l correspond to 
a different situation. All these points intersect the r- 
decaying model tracks at ages > 1 Gyr, so such models 
allow for the underlying old stellar component that is 
present in these galaxies. At the same time, although 
the SFR is exponentially declining, it is still sufficient to 
explain the observed average [Oil] fluxes (cf. Figure fT6]). 
For example, the stacking point for Ltir < 1-7 X 10 11 L Q 
galaxies at D n (4000) = 1.38 and H 5 EW= 3.78 A is con- 
sistent with a r = 0.2 Gyr model at an age ~ 1.4 Gyr. 
By that time, the initial SFR has decayed by a factor 
of e -7 « 10~ 3 . This situation is feasible: a galaxy that 



started its life being a ULIRG would still be producing 
a few M Q /yr 1.4 Gyr later. 

Also, the stacking point at D n (4000) = 1.63 could be 
well explained by the model of a 6 Gyr-old galaxy with a 
much slower SFR decaying time (see the r = 1 Gyr track 
in Figure[T4"|). In spite of being quite old, the SFR of such 
a galaxy is still non-negligible (e -6 rj 2.5 x 10 -3 times its 
initial value). Thus, this model simultaneously accounts 
for the absorption lines and [Oil] emission observed in 
the average spectrum. 

It is important to remark that the fact that T-decaying 
models are able to explain the spectral characteristics 
of these less luminous IR galaxies, the secondary-burst 
scenario cannot be excluded. In all cases, one can find 
the right combination of burst amplitude and onset age 
(i.e. the age of the galaxy when the burst is produced), 
in order to reproduce the average location of galaxies 
with different D n (4000) values in the H 5 EW - D n (4000) 
diagram. 

The IR galaxies for which r-decaying models offer a 
good explanation for the on-going star formation (i.e. 
those with un-corrected D n (4000) > 1.2) constitute ~ 
40% of our Ltir < 3 x 10 11 L Q sample at 0.6 < z < 1.0. 
This result clearly shows that a quiescent mode of star 
formation can drive the IR emission in a significant frac- 
tion of the less luminous LIRGs and IR normal galaxies. 

5.3. Candidates to young starbursts 

Although the average spectra of most 24 (im galaxies 
at 0.6 < z < 1.0 are characterised by having high-order 
Balmer lines such as H^ in absorption, in some objects 
those lines are present in emission. A significant Hi- 
emission is expected for galaxies at the very early stages 
(~ 10 7 years) of a new burst of star formation, when the 
H,5 absorption associated with intermediate-age stars is 
still not sufficient to compensate the emission produced 
by the youngest ones. 

As we explained in C08 and in Section 15.2.21 of this 
present paper, the H5 EW measurements on the spectra 
of our IR galaxies generally have too low signal-to-noise 
ratios as to safely consider them on an individual basis. 
In spite of this fact, we used our tentative Us EW mea- 
surements and secure D n (4000) values to select galaxies 
that are candidates to be hosting young starbursts. 

Our present sample of 1587 faint 24 fim galaxies at 
0.6 < z < 1.0 contains 31 galaxies with D n (4000) < 1.0 
that have H 5 EW> 5 A and < d(H 5 )/H 5 < 1 (i.e. 
H,5 is securely in emission). These galaxies constitute 
18% of all our faint sources with D n (4000) < 1.0 at 
0.6 < z < 1.0, and around 3% of those that are within 
~ 10 s yr of a new burst of star formation (cf. Figure 
[T4"|) . So, the clear cases of galaxies hosting very young 
starbursts constitute a minor fraction of our IR sample 
(we extracted a similar conclusion for brighter IR sources 
in C08). 

There are two alternative explanations for this situa- 
tion. One possibility is that the IR phase is not equally 
likely within the first ~ 10 8 yr of a new burst, i.e. galax- 
ies containing very young bursts are still too faint to be 
detected at IR wavelengths. The other possibility is that 
the earliest stages of star formation are heavily obscured. 
If only the gas had a very large extinction, the observed 
Us emission would be weak, and the H^ EW> 5 A cut 



The optical spectra of 24 /an galaxies in COSMOS: II 



11 



could be missing young starburts. If both the gas and 
stars were very obscured, the galaxies would look faint 
in the optical wavebands and, in that case, they could 
simply be missed by the zCOSMOS-bright I < 22.5 AB 
mag selection. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we have identified 3244 24 /im -selected 
galaxies with 0.06 < S^^m < 0.50 mJy within the 
zCOSMOS-bright 10k sample, and studied different opti- 
cal spectral properties depending on their rcdshifts. Our 
sample spans e.g. total IR luminosities ~ 3 x 10 10 < 
Ltir < 9 x 10 10 L Q at z = 0.5, and - 9 x 10 10 < L T ir < 
5 x 10 11 Ij0 at z = 1. This means that our sources are 
IR normal galaxies at z < 0.5, and a mixture of IR nor- 
mal galaxies and LIRGs at 0.5 < z < 1.0. This works 
complements our previous study of the zCOSMOS-bright 
spectra of more luminous IR galaxies at similar redshifts 
(C08). The fainter luminosity range studied here is par- 
ticularly important, as it contains the bulk of the IR 
background energy at these redshifts. 

For galaxies at 0.2 < z < 0.3, we took advan- 
tage of the simultaneous presence of the Ha A6563 and 
H/3A4861 lines in the zCOSMOS spectra to study the 
relation between line emission, broad-band colours and 
dust extinction. We also discussed the possibility for the 
IR emission to be associated with dust heated by old 
rather than young stellar populations. We summarise 
our results as follows: 

• 92% of our 24 /im galaxies at 0.2 < z < 0.3 have 
Ha A6563 EW> 5 A, which indicates that the vast ma- 
jority of IR sources at these rcdshifts are associated with 
significant on-going star formation. 

In most of these sources, the relation between the gas 
obscuration, i.e. L(Ha)/L(H/3), and IR/UV attenua- 
tion can be reasonably explained with different reddening 
laws of common use in the literature. However, for up to 
16% of our sample at 0.2 < z < 0.3, the Balmer decre- 
ment is much larger than predicted by any of these laws. 
We believe that this is produced by a combination of 
an inhomogenous dust distribution and the fact that the 
zCOSMOS spectra only receive the light from a 1-arcsec 
width slit centred at the source. The presence of a very 
obscured star-forming region towards the centre can pro- 
duce very high L(Ha)/L(H/3) ratios. We note that this 
effect should be much less important at higher redshifts, 
where the 1-arcsec width covered by the VIMOS slit col- 
lects a more representative fraction of the total galaxy 
light. 

It is necessary to emphasize that, although the broad- 
band SEDs of many low-redshift IR galaxies are dom- 
inated by old stellar populations, the optical spectra 
reveal the presence of on-going star formation in most 
cases. Thus, for these galaxies, one cannot conclude that 
the IR light is produced by dust heated by old stars. 
More likely, the young stellar component is the main en- 
ergy source for the IR emission. Of course, this does not 
exclude the possibility that the old stars have a minor 
contribution to the dust heating. 

• Although there is no one-to-one relation between 
the Ha A6563 and H/3 A4861 emission and the galaxy op- 
tical broad-band colours, around 80% of the objects 
with HaA6563and H/3A4861EW > 5 A are in the blue 
(B — i < 1.3) cloud, and around 70% of those with 



H/3 A4861 EW < 5 A are in the 'green valley' or redder 
(B-i> 1.3). 

• The location in the BPT diagram of our galaxies with 
Ha A6563 EW > 5 A suggests that ~ 22% of them arc 
composite systems. However, in none of these galaxies 
can the presence of an active nucleus be confirmed with 
the X-ray or IR AC data available in the COSMOS field. 
This fact suggests that, although an AGN component 
might be present in these galaxies, in any case is not 
dominant, and that aperture effects could also influence 
the fraction of sources identified as composites. 

• We identified only 5 galaxies in our sample at 0.2 < 
z < 0.3 which have no spectral signs of star formation 
and, besides, have elliptical or lenticular morphologies. 
We believe that these are the only candidates in our sam- 
ple at 0.2 < z < 0.3 for which most of the IR emission 
could genuinely be produced by dust heated by old rather 
than new stellar populations. 

The second part of our analysis has been devoted to 
galaxies at z > 0.5. We used the zCOSMOS spectra 
to determine oxygen abundances, measure the D n (4000) 
parameter and the average H,5 EWs. Our results are the 
following: 

• 79% of our galaxies with H/3 A4861 EW > 5 A at 
0.5 < z < 0.7 have degenerate metallicity values. This is 
a consequence of estimating metallicities based on only 
three emission lines (through the R23 parameter) . When 
we place these galaxies in a stellar mass-metallicity dia- 
gram, we see that the upper-branch solution is compat- 
ible with the mass-metallicity relation in the local Uni- 
verse within the error bars. This implies that our data 
do not allow us to conclude on any evolution of the mass- 
metallicity relation between z ~ and z ~ 0.6. Among 
the 21% of galaxies with non-degenerate metallicity mea- 
surements, we find that more than a half are well below 
the mass-metallicity relation. Over our total IR-galaxy 
population at 0.5 < z < 0.7, they put a secure lower 
limit of 12% to the fraction of outliers to this relation. 
This percentage is quite smaller than that observed for 
LIRGs and ULIRGs by C08. However, the comparison of 
our results with those for optically-selected galaxies with 
similar stellar masses and well-determined metallicities 
suggests that the 12% fraction is indeed a lower limit. 

• In contrast to the results for more luminous IR galax- 
ies, we find that our present sample at 0.6 < z < 1.0 dis- 
plays a wide range of D n (4000) values. Still, - 60 (80)% 
of our sources have un-corrected (corrected) D n (4000) < 
1.2. 

• A combined analysis of our present faint IR sample 
and previous C08 brighter sample at 0.6 < z < 1.0 has al- 
lowed us to investigate the possible modes of star forma- 
tion in IR galaxies. The burst- like mode of star formation 
is necessary to simultaneously explain the different aver- 
age spectral properties of the most luminous IR galaxies 
(Ltir > 3x 10 11 L Q ). Instead, for nearly 40% of the less 
luminous LIRGs and most luminous IR normal galaxies, 
a continuous, quiescent star-formation mode is possible. 
Realistically, the Ltir for the transition at 0.6 < z < 1.0 
should be considered to be (3 ±2) x 10 11 L Q , to take into 
account the errors in the determination of the total IR 
luminosities and the data binning. 

Our findings are consistent with other recent results in 
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the literature, which indicate that around a hall ol IR 
galaxies ol intermediate luminosities at z ~ 1 are typical 
disks that do not show any sign ol disturbance (see e.g. 
Elbaz et al. 2007). They could also be explained in the 
scenario proposed by theoretical studies that suggest that 
some galaxies could have had a star formation history 
characterised by steady gas flows since early times (Dekel 
et al. 2009). 

In summary, our optical spectral analysis has clearly 
shown the diverse nature ol the galaxies that make the 
bulk ol the mid-IR background at z < 1. And, even 
when star formation seems to be the most important 
driver ol the IR emission at these redshilts, IR galaxies 
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Fig. 1.— The zCOSMOS-bright 10k and the SCOSMOS/MIPS shallow and deep sample coverage fields. 
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Fig. 2. — The zCOSMOS redshift distributions of the new 24 fim sources identified in this work, and those considered for spectroscopic 
analysis (empty and filled histograms, respectively). 
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Fig. 3. — HoA6563and H/3A4861EWs versus luminosities for our IR, galaxies at 0.2 < z < 0.3 (left and right-hand panels, respectively). 
All measurements are corrected for stellar absorption, and line luminosities also include aperture corrections. 
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Fig. 4. — Left: the distribution of observed (B — i) colours for 24 fim -detected and all zCOSMOS galaxies at 0.2 < z < 0.3. Right: the 
(B — i) colour distributions for the 24 /an galaxies segregated by their Balmer-line EWs. 
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Fig. 5. — The observed (B — i) colour versus the Balmer decrement for our 24 fim galaxies with HqA6563EW >5Aat0.2<z< 0.3. 
Typical error bars for galaxies with H/3 A4861 EW > and < 5 A are shown at the top of the plot (with a circle and a diamond, respectively). 
The vertical dashed line indicates the Ho A6563 /H/3 A4861 ratio expected in a case B recombination with T = 10, 000 K (Osterbrock 1989). 
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Fig. 6. — The average stacked spectra of 0.2 < z < 0.3 24 «m galaxies with individual HaA6563EW > 5 A and H/3A4861EW < 5 A, 
according to their (B — i) colours. 
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Fig. 7. — The average stacked spectrum of 24 /im galaxies at 0.2 < z < 0.3 with individual HaA6563EW < 5 A. 
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Fig. 8. — Broad-band SEDs of galaxy candidates for having dust heated by old rather than young stars (0.2 < z < 0.3). Circles indicate 
rest-frame total broad-band fluxes (note that the error bars are smaller than the symbol size). Lines of different style correspond to 
renormalised SED model templates from Polletta et al. (2007), constructed using the GRASIL code (Silva et al. f998) and empirical tR 
spectra: Elliptical (solid), SO (dashed) and Sc (dot-dashed). 
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Fig. 9. — The Balmer decrement versus the far-IR/UV attenuation. Symbols are explained in the plot label. Right-pointing arrows 
indicate lower limits for objects non-detected in t he N UV. Plus-like symbols within the circles/diamonds indicate objects that are in the 
composite region of the BPT diagram (see Figure [TO) . Typical error bars are shown on the left-hand side of the plot. Lines of different 
styles indicate the relation predicted with different reddening laws: SMC (Prevot et al. 1984; double-dot-dashed); MW (Fitzpatrick 1999; 
dashed); Calzetti et al.( 2000; solid) and a composite law (dot-dashed). The horizontal dashed line shows the Ha A6563 /H/3 A4861 ratio 
expected in a case B recombination with T = 10, 000 K (Osterbrock 1989). 
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Fig. 10. — The BPT diagram for our 24 (im sources at 0.2 < z < 0.3. Symbols are the same as in figure l9l The dashed line indicates the 
starburst/AGN division line proposed by Kewley et al. (2001), while the dot-dashed line corresponds to the limit for composite systems 
(Kauffmann et al. 2003a). 
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Fig. 11. — The [OIII]/[OII] ratio versus the R23 parameter for the galaxies in our metallicity sample at 0.5 < z < 0.7. Small and large 
circles refer to galaxies that have degenerate and non-degenerate metallicity estimates, respectively. 




Fig. 12. — Oxygen abundances versus assembled stellar masses for our 24 /im galaxies with H/3A4861EW > 5 A at 0.5 < z < 0.7 (filled 
circles: large and small for non-degenerate and degenerate O/H values, respectively). The thick blue crosses indicate the lower-metallicity 
alternatives of the blue circles. The plus-like symbols inside the circles refer to X-ray-detected AGN. The average error bars shown in 
opposite corners of the plot correspond, respectively, to metallicities in the upper and lower branches (non-degenerate metallicities around 
12 + log 10 (O/H) ~ 8.6 have similar errors to the lower-branch metallicities). Open circles and thin black crosses correspond to the more 
luminous 24 (im galaxies analysed in C08, while open triangles indicate CFRS galaxies with known metallicities (Maier et al. 2005). The 
solid line shows the best-fit mass-metallicity relation obtained by Kewley & Ellison (2008) on SDSS galaxies at 2 < 0.1. The solar 
metallicity derived by Asplund et al. (2004) is also shown for reference (dotted line). 
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Fig. 13. — Rest-frame 24 fim luminosities versus D n (4000) parameter values for our 24 (im galaxies at 0.6 < z < 1.0. Filled circles 
correspond to galaxies in our current faint 24 fim sample, while open circles indicate the C08 24 (i m so urces. The dashed lines delimit the 
bins we have adopted for spectral stacking and H<5 EW versus D n (4000) measurements (cf. Figure [T4t . The D n (4000) values in this figure 
are not corrected for extinction. 
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Fig. 14. — The H<5 EW versus D n (4000) diagram. The average location of our zCOSMOS 24 /im galaxies at 0.6 < z < 1.0, classified in 

intervals of mid-IR luminosity ^L^ 4Mm and D n (4000) values, are indicated by filled squares (note: the average D n (4000) values shown in 
this plot are statistically corrected for dust extinction; see text). The total IR luminosities Z/tir shown in the labels have been obtained 
applying the conversion formula Ltir = 6856 X {uL^ ' Jm )°' 71 (Bavouzet et al. 2008). The thick lines show the tracks produced by various 
exponentially-declining star formation histories: from top to bottom, r = 0.01, 0.2, 1, 2Gyr and r — » oo (i.e. constant SFR). Tracks only 
show the evolution of each model up to 6 Gyr, which is approximately the age of the Universe at z = 1. The filled circle intersecting each 
track marks the position at 1 Gyr. Thin solid lines indicate the effect of secondary bursts of different amplitudes (the relative amplitude 
with respect to the primary burst is indicated with a percentage), produced at different ages of an underlying old stellar population. 
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Fig. 15. — Stacked spectra for the less luminous IR galaxies in our sample, and the smallest D n (4000) values (the bins displayed in the 
labels are those used for the stacks, i.e. D n (4000) is not corrected for extinction). 
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Fig. 16. — The average [Oil] fluxes versus D n (4000) values measured on our stacked spectra at 0.6 < z < 1.0, corresponding to different 
i/Ljr' im and D n (4000) bins. Symbols are the same as in Figure [141 



